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ABSTRACT 

We provide predictions on small-scale cosmological density power spectrum from supernova 
lensing dispersion. Parameterizing the primordial power spectrum with running a and running 
of running /3 of the spectral index, we exclude large positive a and /3 parameters which induce 
too large lensing dispersions over current observational upper bound. We ran cosmological N- 
body simulations of collisionless dark matter particles to investigate non-linear evolution of 
the primordial power spectrum with positive running parameters. The initial small-scale en¬ 
hancement of the power spectrum is largely erased when entering into the non-linear regime. 
For example, even if the linear power spectrum at fc > lOhMpc”^ is enhanced by 1 — 2 orders 
of magnitude, the enhancement much decreases to a factor of 2 — 3 at late time (z ^ 1.5). 
Therefore, the lensing dispersion induced by the dark matter fluctuations weakly constrains 
the running parameters. When including baryon-cooling effects (which strongly enhance the 
small-scale clustering), the constraint is comparable or tighter than the PLANCK constraint, 
depending on the UV cut-off. Further investigations of the non-linear matter spectrum with 
baryonic processes is needed to reach a firm constraint. 

Key words: cosmology: theory - gravitational lensing: weak - large-scale structure of the 
universe - cosmological parameters - inflation 


1 INTRODUCTION 


The plethora of cosmological observations has turned Cosmology 
into a quantitative science. From combining several probes that ob¬ 
serve the Universe at different epochs and have different systemat- 
ics and statistics, emerged the ‘Concordance Model’ of Cosmology, 
a six parameter model, most of them measured to the accuracy of a 
percent. Among these probes, type la supernovae (SNe) are a pow¬ 
erful cosmological tool to directly measure the expansion histor y 
of the Universe (see e.g. a recent review of lWeinberg et al.ll2013h . 
The type la SNe are known as standard candles and one can mea¬ 
sure the luminosity distances to them accurately. Several SNe sur¬ 
veys have set strong constraint s on cosmological models from the 
distance-redshift relation (e.g. ^esset^3li22i Pgdmutteret^^ 


ll999l:lRiess et al.l2007l:ISullivan et aLlioIIl: ICampbell et al .120131) . 

Flowever the measured apparent magnitudes, have a residual scat¬ 
ter arising from its intrinsic scatter and effects due to line-of-sight 
(LOS) structures. The intrinsic scatter (~ 0.4 mag) can be signif¬ 
icantly reduced to ~ 0.1 mag by empirically calibrating the lumi¬ 
nosity curves. The scatter due to photon deflection along the line of 
sight is composed of many different physical effects. The dominant 
effects are peculiar velocities at 2 < 0.1 and gravitational lensing at 
2 > 0.3. Within the realm of cosmological perturbation theory and 
the stochastic nature of the LOS structures, all effects are expressed 
as an integral over the power spectrum with appropriate kinematical 


factors teen-Davan et alj|2012h . The lensing effects create residu¬ 
als from the best-fitting curve in the magnitude-redshift relation. 
The lensing magnifications of SNe can be extracted by correlating 
the residuals with the surface densities of nearby foreground galax¬ 
ies. The lensing dispersion is roughly proportional to the SNe red- 
shift 2 as ~ 0.002 mag (e.e. lHolz & Lindej2005l : lBen-Davan et alj 
l2013h . 


Detection of the lensing magnifications has been reported by 
several groups jjonsson et al.ll200^ l20ld : iKronborg et al.llToioi : 
iKarnenka et alj2o'l3 : Smith et alj2014 ). Jonsson et al. (2010, here¬ 
after JIO) used a sample of 175 SNe at 0.2 < 2 < 1 from the 
SNLS 3yr data, and measured SNe lensing magnifications by corre¬ 
lating the SNe fluxes and its foreground nearby galaxies (modeled 
as singular isothermal spheres). Modelling the lensing dispersion as 
Cm(zs) = Bzs, they obtained a constraint B ~ O.OSSIq q 4 ® mag 
(la) and an upper bound B < 0.095(0.12) mag with l(2)-sigma 
confidence. More r ecently the Joint Lightcurve Analysis (JLA) 
teetoule et al.ll201^ combined SNe from various observations to 
compile 740 SNe for cosmological parameter inferenc e, yielding 
fJm = 0.295 ± 0.034 very close to the PLANCK result jAde et all 
l2015h . The JLA did the following errors analysis. First it consid¬ 
ered the peculiar velocity dispersion, that is relevant only for very 
low redshift. In addition, it used JIO central value for the lensing 
dispersion. Last, it considered all additional possible sources of er¬ 
ror and named it ‘coherent dispersion’, acoh- Adding all of these 
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in quadrature produced the total dispersion, atotiz 1) ^ 0.12 
in magnitude. There is also a clear trend of the total dispersion de¬ 
creasing with redshift, Fig. 7 in JLA, suggesting that lensing might 
even be smaller than JIO central value. □ Throughout the paper we 
will take this number ^ 0.12 as a conservative upper bound 
of the lensing dispersion, and ^ 0.095 as a ‘plausible’ upper 
bound, based on Icr of JIO. 

The photons’ scatter biases the measurement. However, given 
a cosmological model, we can use gravitational lensing as a probe 
of Cosmology and Astrophysics. There are several advantages of 
using SNe lensing as a cosmological probe. First, it has differ¬ 
ent systematics than other probes. Second, it is sensitive to dif¬ 
ferent scales. For example, we shall see that lensing dispersion is 
mainly sensitive to 10“^/iMpc“^ ^ ^ 10®/iMpc“^. Third, 

even current data (or upper bound) is sufficient to constrain the 
cosmological parameters. The gravitational lensing of type la SNe 


has been suggeste d as a probe of cosmological abundance of com- 


pact objects (e.g. RaucfJ 

1991: Holz 1998112001: Metcalf & SilkI 

19991: Mortsell et alJ2001 

: Mintv et al J2002l; Metcalf & Silkl2007l; 

Yoo et alj|2008h and for measuring the present amplitude of mat- 


ter densitv fluctuations (e.g. Metcalf] 19991: Hamana & Futamase 

200C 

: Valagead 20001; Dodelson & Vallinotto 

20061; Marra et al 

2013 

:IOuartin et alJ 

2014; Fedeli & Moscardini 

l20l4llwata & Yoc 

2015 

). For instance. 

Castro & OuartinI J2014h recentlv measured as 


(the present amplitude of density fluctuations at 8/i“^Mpc) using 
the non-Gaussianity of SNe magnitude distribution from 700 SNe 
at 2 < 0.9. 


Recently. iBen-Davan & Kalavdzhv^ l l2014h used the lensing 
dispersion to constrain small-scale amplitude of primordial cur¬ 
vature perturbation. They modelled the primordial power spec¬ 
trum using running a = dris/dlnk and running of running 
/3 = d^Us/dlnk^ of the spectral index, and then discussed an 
observational constraint on a and j3 from the lensing dispersion. 
A subsequent work consi dered other param eterizations of the pri- 
mordial power sp ectrum, ( lBen-Davanll2014h . Current observations 
jAde et alj|20ll) measure ns{ko) — 0.96, for the parameteriza¬ 
tion Ph = where ko is a suitable “pivot scale”. 

Focusing again on “running”, we consider: 


ns{k) 


ns{ko) + 


ajko) Pjko) 

2 ko^ 6 



( 1 ) 


There is no compelling evidence for the more generalized form, but 
the best Ax^, that takes into account the suppression of low multi¬ 
ples in the temperature anisotropies spectrum is with a = 0, /? = 
0.029 jAde et alJboTsh FI Given a single canonical scalar field 0 
and a potential V{(/)), the above observables can be described by 
the slow-roll parameters, or the potential and its derivatives at the 


^ While the lensing dispersion grows with redshift, it is by no means the 
only contribution to the dispersion that can evolve with redshift. For exam¬ 
ple the SNe population can in principle evolve with redshift degrading the 
efficiency of SNe lensing as a probe. We thank Adam Riess for raising this 
issue. 

^ Non-trivial running i s a p rediction if e is non monotonic, see 
Ben-Dayan & Brustei j Eoiol). and a Higg s infl ation version in 
Hamada et al H20l4) : lBezrukov & Shaposhniko^ j20l4) . 


’CMB point’: 

Us = 1 + 2ri — 6e ( 2 ) 

Q = 1667? - 24e^ - 2f (3) 

P = -192e^ -f 192dri - 32er?^ - 24e^^ -f 2r]f + 2u;® (4) 

1 /V'V V" 2 V'V" 3 

where prime denotes derivative with respect to the inflaton. The 
best constraints on a, /3 with fcp = 0.05Mpc~^, n^ffco) ~ 
0.96 are given by PLANCK jPlanck Collaboration et alj [201461 : 
lAde et all l2015h . These analyses are only probing the range 
Ho ^ k < lMpc“^. The lensing dispersion, cTm is sensitive to 
0.01Mpc“^ < k < 10® Mpc“^, thus giving access to 2 — 3 more 
decades of the spectrum. Hence, amiz) is particularly sensitive to 
the quasi-linear and non-linear part of the spectrum. In terms of in¬ 
flation, the direct measurement of fc < lMpc“^ corresponds to 
about 8 e-folds of inflation, leaving most of the power spectrum of 
~ 60 e-folds out of reach. Hence, even after PLANCK there is still 
an enormous space of inflationary models allowed. It is therefore of 
crucial importance to infer as much of the spectrum as possible for 
a better inflationary model selection. The lensing dispersion con¬ 
strains additional 4 — 7 e-folds, yielding a total of 12 — 15 e-folds. 

The main limit a tion of the a n alysis in 
iBen-Davan & KalavdzhvanI j2014l) and iBen-Dav^ j2014l) . was 
the knowledge of the non-linear power spectrum k > l/iMpc“^ 
at late times 2 < 2, in the case of running a and running of 
running /3. In the absence of numerical simulations with proper 
initial conditions, enhancement of power on non-linear scales was 
parameterized by the enhancement that occurs in the Concordance 
Model (i.e. a = /3 = 0). In this paper, we overcome the main 
limitation of previous works. We ran N-body simulations with 
non-zero a, P and derive a modified HaloFit formula for the 
non-linear power spectrum. We calculate the lensing dispersion 
and compare it with the observational upper bound atot ^ 0.12. 
Then, large positive running parameters a and P can be rejected. 
We will also include the baryonic cooling effects, which strongly 
affect the small-scale clustering, in our investigation. 

Before outlining the rest of the manuscript a valid question 
arises regarding the logic of considering just the lensing dispersion, 
rather than the full lensing probability distribution function (pdf). 
The calculated and experimentally quoted lensing dispersion is that 
of the distance modulus, m ~ —2.5 log^Q F . The flux F is usuall y 
simplified to have a log-normal distribution jHolz & Lindeill200^ . 
In that case, the distance modulus is a Gaussian random variable 
and all the information is contained in the dispersion, and one does 
not need to consider higher moments of the pdf. Consideration of 
higher moments is relevant for deviations from the log-normal dis¬ 
tribution of the flux and is beyond he scope of this work. 

The paper is organized as follows. In sec.2, we briefly describe 
the lensing dispersion formula. In sec.3, we first present the primor¬ 
dial power spectrum parameterized with the running parameters (a 
and P), and then run the cosmological N-body simulations to fol¬ 
low its non-linear evolution accurately. The fitting formula of the 
non-linear matter power spectra are presented in this section and in 
the Appendix. Sec.4 presents our main results of the lensing dis¬ 
persion and observational constraints on the running parameters. 
Here, we also include the baryon cooling effects. Conclusions and 
a discussion are in section 5. Throughout this paper we adopt the 
concord ance ACDM model, which is consi stent with the PLANCK 
results jPlanck Collaboration et al.ll2014ah . The model is charac¬ 
terized by the matter density Cm = 0.3134, the baryon density 
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f2b = 0.0486, the cosmological constant density f^A = 0.6866, 
the amplitude and the spectral index of primordial curvature power 
spectrum As = 2.20xl0~®, Ua = 0.96, and the Hubble expansion 
rate today Hq = 67.3 km s“^ Mpc“^. We have verified that vary¬ 
ing the other cosmological parameters does not change the conclu¬ 
sions beyond the quoted error bars (see lSen-Davan & Kalavdzhvarj 

[20Ti) . 


2 LENSING DISPERSION 


When a distant SNe is gravitationally lensed by intervening struc¬ 
tures, its flux is multiplied by a magnification factor /i and its 
apparent magnitude (or distance modulus) is changed by 5m = 
—2.5 logj^Q /i. For a perturbed FLRW Universe, in the Poisson 
gauge, one starts with the line of sight (LOS) first order lensing 
contribution to the distance modulus. 


5m{rii°f 


5 drii rji - rii°^ ^ ^ 
In 10 7 ^( 0 ) Ar] r]o -r/i ^ 


(5) 


where the gravitational potential is evalu¬ 

ated along the past light-cone at Vi = rjo — rji, rjo is the 
observer conformal time, is the conformal time of the 

source with unperturbed geometry, Arj[z) = rjo{z) — ff\z) 
= dji / -f v)^ + Ha and A 2 is the 2D a ngul^ 

Laplacian, see Isen-Davan et alj l l2012h : ISen-Davan et alj l l2013h 
for technical terms and further explanations. Squaring © and tak¬ 
ing the ensemble average in Fourier space at a fixed observed 
redshift, gives the variance, (j^{z). The double LOS integral at 
2 > 0.3, is dominated by the equal time part, and with Si(a: 

1) « 7 r/ 2 , gives: 

j dkPflk,nCe ( 6 ) 



Figure 1. Dimensionless curvature power spectrum 'p 7 ^(fc) in our running 
spectral index models: (o,/?) = (0.06,0.06), (0,0.06), (0.03,0.04), 
(0.02, 0.02), (0, 0.02), (0.02, 0), (0, 0) from top to bottom. Here, in ver¬ 
tical axis, the power spectrum is divided by its amplitude As. 


iMortsell et al.l2Q0lh to the present Horizon scale Hq ^. From equa¬ 
tion 0, it is clear that large scales k < fceq have a negligible con¬ 
tribution. At smaller scales, k ~ l/iMpc“^, perturbation theory 
breaks down. Additionally, we expect non-Gaussian effects to be¬ 
come relevant. Last, the Nyquist wave number of the simulations 
soon to be discussed is at most fcNyq = 64.3/iMpc“^. We therefore 
limit the integration at most to = lO^/iMpc”^, while fcmin = Hq. 

In the following sections, we will evaluate the non-linear mat¬ 
ter power spectrum using cosmological N-body simulations. Here, 
the (dimensional) matter power spectrum Pm and the (dimension¬ 
less) gravitational potential power spectrum P-^ are related via, 


where P.i< is the dimensionless power spectrum of the gravitational 
potential. Hence, the lensing dispersion of supemovae is a direct 
measurement of the integrated late-ti me power spectrum. Le t us 
stress that it is both UV and IR finite teen-Davan et alj|2013h . In 
general, the k^ enhancement makes cTm a sensitive probe to the 
small scales of the power spectrum. This sensitivity is manifest if 
we switch to dimensionless variables, t) = Hop and p = kjkao- 


2 

fTm — 


In 10 


TT 

Aip 


J dpi j dpP^(p,pi)p^ (7) 

x(»7i - nff^ipo - pif ■ 


It is easy to see that: (a) The relevant physical scales are Ho and 
fceq, which give an enhancement of (fceq/Po)^- ( 6 ) The disper¬ 
sion is really sensitive to scales smaller than the equality scale 
p > 1. (c) The nonlinear power spectrum (NLPS) has an addi¬ 
tional redshift dependent physical scale which is the onset of non¬ 
linearity A:nl. For a given redshift, parameterizing the NLPS as 
~ C{k/k^C)'^, V < —3, from some /cnl, will have an additional 
parametric enhancement of (feNu/feeq)^ to the lensing dispersion. 
Note that (Jm is in units of magnitude and the dispersion in equa¬ 
tions 0 and m. is exactly the standard resu lt in weak gravitational 
lensing (e.g. lSartelmann & Schneide3l200lh . In principle, the lens¬ 
ing dispersion can probe the density fluctuati ons from a size of SNe 
at its maximum luminosity (« 10^® cm) (e.g. lMetcalf & SilkI 19991 : 


P^{k,z) = .^Htnl,{l + zfk-^Pm{k,z). ( 8 ) 


3 NON-LINEAR MATTER POWER SPECTRUM 

We run cosmological N-body simulations to investigate the non¬ 
linear matter power spectra of the various running spectral index 
models. Our purpose is to obtain a fitting formula of the non-linear 
power spectra and calculate the lensing dispersion in equation 0 . 


3.1 Linear power spectra with running spectral Indices 

The dimensionless curvature power spectrum is given by 

P^W-Aa(-) . (9) 

where ko = 0.05Mpc“^ is the PLANCK pivot scale, As = 
2.20 X 10“® is the amplitude and ris = 0.96 is the spectral index. 
For the running a and the running of running /3 of the spectral in¬ 
dex, we adopt seven models as (a, P) = (0, 0), (0, 0 . 02 ), ( 0 . 02 , 0), 
(0.02, 0.02), (0.03, 0.04), (0, 0.06) and (0.06, 0.06). Note that As 
is fixed for the running models (a 7 ^ 0, /? 7 ^ 0). Then, the present 
mass variance as at 8 / 1 “ ^Mpc weakly increases as the running pa¬ 
rameters (a and /?) increase: as = 0.829 for (a, P) — (0,0), 
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0.831 for (0,0.02), 0.833 for (0.02,0), 0.835 for (0.02,0.02), 
0.839 for (0.03,0.04), 0.834 for (0, 0.06), 0.847 for (0.06, 0.06), 
respectively. FigUplots the curvature power spectrum in equation 
© in our running models. The small-scale enhancement of Vnik) 
is clearly prominent. Large a,/3 models strongly enhance Vnik) 
by an order of magnitude for fc > 10 /iMpc“^. The PLANCK re¬ 
sults already excluded large ru nning models of (a, 0) — (0, 0.06 ) 
and (0.06, 0.06) at 2a level jPlanck Collaboration et aljr2014bh . 
However, it is useful to study the general feature of the non¬ 
linear gravitational evolution with such strong small-scale enhance¬ 
ment. Furthermore, it provides an independent crosscheck, to the 
PLANCK analysis, and as we have mentioned, it probes scales be¬ 
yond PLANCK’S reach. 

We prepare the linear matter power spectrum with the running 
parameters, Pm,iin(fe; ce, P), as follows. We first evaluate the linear 
power spect rum without the ru nning Pm,iin(fc; a = 0, /3 = 0) us¬ 
ing CAMS jLewis et al.ll2000h . Then, simply multiplying a factor 
from equation (|9) to the output of CAMB as, 

^’m,lm(fc; a, P) = f’m,lin(fc; Q = 0, /3 = 0) 

1 . \ (a/2) ln(fc/feo)-|-(/3/6)[ln(fe/fco)]^ 

■ ( 10 ) 



3.2 Cosmological N-body simulations 

We run cosmological numerical simulations of dark matter parti¬ 
cles in cubic boxes with a side of 100/i“^Mpc. In order to check 
the numerical resolution, we run high- and low-resolution simula¬ 
tions in which the number of particles are 2048^ and 1280^. The 
initial matter power spectrum is prepared following the procedure 
in section [TT] The initial positions and velocities of particles are 
given at redshift gjnit = 99 based on second-order Lagrangi an per¬ 
turbation theory dCrocce et alJl200^ : lNishimichi et alj2009h . 

To follow the gravitational evolution of the dark matter 
particles, we employ a pu blicly available tree-PM code GreeM 
dishivama et al.ll20(Ml2012h . GreeM is tuned to accelerate the tree 
gravitational calculation, and it is fast especially in the strongly 
non-linear regime. The PM meshes are set to be 1024^(640®) and 
the particle Nyquist wavenumbers (which is determined by the 
mean particle separation) are fcNyq = 64.3(40.2)/iMpc“^ for the 
high-(low-)resolution simulations. The small-scale enhancement 
seen in Fig[T]is included up to the Nyquist wavenumber (k < AiNyq) 
in the initial conditions of our simulation. In other words, the ini¬ 
tial linear Piim(fc) is zero at A: > ZtNyq. The gravitational softening 
length is set to be 3% of the mean particle separation. The sim¬ 
ulation snapshots are dumped at redshifts z — 0, 0.3, 0.6, 1 and 
1.5. We prepare 3(6) independent realizations for the high-(low- 
)resolution simulations for each (a, P) model to reduce the sample 
variance. Our simulation settings are summarized in Table[T] 

Fig[2]shows our simulation results of the matter power spectra 
with the various running parameters (a, p) atz = 0, 0.3, 0.6, 1 and 
1.5. The blue diamonds (black filled circles) with error bars are the 
mean power spectra with the standard deviations obtained from the 
realizations of high-(low-)resolution simulations. The blue (black) 
dotted lines are the shot noise for the high-(low-)resolution simu¬ 
lations. The high-resolution results are slightly higher (about 10%) 
than the low-resolution ones especially at low redshifts a ^ 0.6. 
The initial Pm,iin(A:) is exponentially enhanced for some large a 
and P, and hence the simulation results would strongly depend 
on the resolution especially at wavenumbers comparable to or 
larger than the Nyquist wavenumber. However the differences be¬ 
tween the high- and low-resolution results are smaller than 20% 


for fc = 3 — 300fcMpc“^ at z = 0 — 1.5. The red solid curves 
are our fitting formula based on the Halofit mo del for the ACDM 
model dSmith et al.ll2003l : lTakahashi et al.ll2012h . but slightly mod¬ 
ified for the running spectral index models. The Halofit agrees 
with the simulation results within 18.3(28.3)% relative error for 
fc < 300(600)/iMpc“^ at 2 = 0 — 1.5. Details of the fitting pa¬ 
rameters are given in Appendix A. 

As seen in the figure, the matter power spectrum Pm(fc) with 
larger running parameters is further enhanced at small scales. For 
example, Pm(fc) with {a, P) = (0.06,0.06) and (0.03,0.04) 
are about 3 and 2 times larger than that without the running at 
fc = 50/iMpc“^. However, the enhancement is much less promi¬ 
nent than those seen in the linear power spectrum (see Fig[Tll 
which shows 1 — 2 orders of magnitude enhancement at fc > 
lO/iMpc”^. Our numerical results show that the small-scale en¬ 
hancement in the initial Pm(fc) are largely smeared out at late 
times. This result exhibits power transfers from large to small 
scales, via mode coupling, and thus suggests that the initial large- 
scale power mainly determines the non-linear small-scale power 
at late times. These power transfers in the non-linear gravita¬ 
tional evolution has been well studied for the models in which 
an ini tial Pni(k) has a peak ( like a Gaussian) at a specific scale 
(e.g., Bagla & Padmanabha nI 19971: Padmanabhan & Ray 200^ 


Bagla & PrasadI l2009l : iNevrinck & Yand l2013l : iNishimichi et alJ 
2014h and in models where the initial Pm(fc) i s damped at small 


scale such as warm dark matter models (e.g., [Little_et_alJ 


White & Cro^ 120001: ISmith & Markovid 1201 iL IViel etahl 


1991 


2012 


Inoue et alT 2015h .~Such a smearing of small-scale enhancement 


makes it difficult to discriminate the running parameters (a, P) 
from the nonlinear Pm(fc). 

Before closing this section, let us comment on effects of 
the primordial small-scale enhancement on a halo mass function. 
We compute the mass function by identifying haloes in the low- 
resolution simulations using standard friends-of-friends algorithm 
with linking length of 0.17 times the mean particle separation. 
Then, as expected, the strong running models predict more numer¬ 
ous small halos, but the differences of the mass functions among 
the seven models are only less than a factor 2 at 2 = 0 — 1.5. 
For example, even for the smallest halo of 4 x 10^h~^M q (which 
contains 100 simulation particles), the largest running model (a = 
P = 0.06) has about 2(1.4) times more haloes comparing to the no 
running model (a = /I = 0) at 2 = 1.5(0). This result suggests 
that the primordial small-scale enhancement does not significantly 
affect the halo mass function at late time (z ^ 1.5). 


4 RESULTS 

This section presents our main results of an observational upper 
bound on the running parameters (a, p) from the leasing disper¬ 
sion. First, subsection l4. II gives the results based on the non-linear 
matter power spectra of dark matter fluctuations obtained in previ¬ 
ous section. Then, subsection 14.21 includes baryonic effects on the 
matter power spectra and the leasing dispersion. We will show that 
the baryonic effects strongly enhance the leasing dispersion and 
tighten the upper bound on a and p. 


4.1 Dark Matter Only 

Fig[^shows our results of the leasing dispersion for the various run¬ 
ning parameters {a,P). Here, we numerically calculate the leas¬ 
ing dispersion in equation © using our fitting formula of Pm(fc). 
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(a,|3)=(0,0) 



(a,P)=(0.02,0) 





T 




{a,P)=(0.02,0.02) 



♦ high-resolution sim. 

• low-resolution sim. 

- modified haiofit model 

z=0, 0.3, 0.6, 1, 1.5 
from top to bottom 

gray symbols are results 
for (I'it) tO.fii 


10 ° 10 ^ 10 ^ 10 ° 10 ° 10 ^ 10 ° 10 ° 10 ° 10 ' 10 ° 10 ° 


k (h/Mpc) 


Figure 2. Non-linear matter power spectra with running a and running of running /9 of the spectral index. The seven panels con'espond to [a, ft) = (0, 0), 
(0, 0.02), (0.02, 0), (0.02, 0.02), (0.03, 0.04), (0, 0.06) and (0.06, 0.06), respectively. The redshifts are z = 0, 0.3, 0.6,1 and 1.5 from top to bottom in 
each panel. In the vertical axis, the power spectrum is multiplied by a factor to show clearly our results. The blue diamonds (black filled circles) are the 
high-(low-)resolution simulation results. The dark gray symbols are the results for the no-running model (a, 0) = (0, 0). The blue (black) dotted lines denote 
the shot noise in the high-(low-)resolution simulations. The red curves are our modified haiofit model (see main text and Appendix A for detailed discussion). 


Table 1. Our Simulation Setting 



L(/i-iMpc) 

/V3 

fcNyq(hMpC“l) 

mp{h~^ Mq) 

z 

N, 

high-resolution 

100 

2048® 

64.3 

1.0 X 10'^ 

0,0.3, 0.6,1,1.5 

3 

low-resolution 

100 

1280^ 

40.2 

4.1 X 10'^ 

0,0.3, 0.6,1,1.5 

6 


The different columns are the side length of simulation box L, the number of dark matter particles Nf the Nyquist wavenumber IcNyq, the particle mass rrip, 

the redshifts of the simulation outputs z, and the number of realizations . 


The figure clearly shows the lensing dispersion strongly depends 
on the running parameters. For instance, the lensing dispersion for 
the strongest running model (a = 0 = 0.06) is about 40% higher 
than the Concordance Model (a = /3 = 0) at 0 < 2 s ^ 2. The bot¬ 
tom panel is the lensing dispersion divided by the source redshift, 
which shows amiza) is roughly proportional to Zs at Za > 0.8 (e.g. 
IfIoIz & Lindeill200a : lBen-Davan et alj2013l) . 

We comment on the interval of the integration of the lensing 
dispersion. The integrand is proportional to Pjn{k)k^ per dink 
from equations and and thus the contribution to the in¬ 
tegral comes from the peak of which corresponds to 

k ~ lO/iMpc”^ from Fig[^ Thus, the lensing dispersion does 
not strongly depend on the interval of the integral, if fcinin(A:max) 
is much smaller (larger) than « 10/iMpc“^. For example, if we 


set the maximum wavenumber fcmax = instead 

of lO^/iMpc”^, the lensing dispersion becomes 6% larger for 
a = P = 0.06 and less than 3% larger for the other models at 
Zs = 1. 

Fig|4]is a contour plot of the lensing dispersion at 2 s = 1 
in a/3 plane. Flere we apply our fitting formula of Pm(fc) in a 
range of 0 < a,/3 < 0.1. The thin black lines correspond to 
(Tm = 0.07, 0.08, 0.09, • • • (mag). The thick line (CTm = 0.095) 
corresponds to the observational upper bound (68% C.L.) from 
the Supernova Legacy Survey (SNLS, JIO). White contours are 
P LANCK constraints of 1 and 2 sigm a regions (Planck-l-WP-l-BAO 
in IPlanck Collaboration et alj|2014d) . The dependence of a,/3 in 
cTm (seen in the solid lines) comes from the fitting formula ( lA4b . 
3.32a + 4.72/3 = const. Then, the SNLS constraint can be rewrit- 
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Figure 3. Lensing dispersion as a function of source redshift for the 
vaiious running a and running of running (3 of the spectral index. Top panel 
is the lensing dispersion in units of magnitude, and bottom one is the lensing 
dispersion divided by the source redshift. 


Figure 4. Contour plot of the lensing dispersion at 2s = 1 in a/3 plane. 
The black solid lines con'espond to crm = 0.07, 0.08, 0.09, • • • in steps 
of 0.01 (mag). The thick line (cjm = 0.095) denotes the observational 
upper bound (1 sigma) from the Supernova Legacy Survey (JIO ). The white 
contours are the PLANCK jPlanck Collaboration et al.ll2014bl) constraints 
of 1 and 2 sigma regions. 


ten as a + 1.42^ < 0.227 (95% CL), and a + 1.42/3 < 0.167 
for the ’plausible upper bound’ ( 68 % CL). Our constraint may 
seem weaker than the PLANCK results, but it is important to 
note, that unlike CMB results, ours extend all the way to /Cmax = 
10^/iMpc“^. So while the result is not compelling as one may have 
hoped for, it still contains new information about the small scales 
of the primordial power spectrum in that it extends PLANCK con¬ 
straints to the /cmax scale, if lensing is actually observed. 

Our numerical results of the lensing dispersion at 0.6 ^ 2:3 ^ 
2 and 0 ^ a, ^ 0.2 can be fitted as, 

- = h{zs) + fi{z^) {x 3-f 2 {zs)x ] mag 

Zs 

( 11 ) 

where x = 3.32a + 4.72/1 and 

fo{zd = 0.061 + 0.010(^5 - 1 ) - 0 . 012 ( 2 s - l)^ 

/i(2s) = 0.038 + 0.003(2s-1), 

f2{zs) = 0.75 + 0.37(2s - 1) - 0.20(2s - if. (12) 

The fitting function can reproduce the our numerical results within 
a relative error of 4.3% for ffm/zs < 0.3mag. Since am is roughly 
proportional to 2 a, 2 s-dependencies of the functions /o,i, 2 (zs) are 
relatively weak. 



Figure 5. The baryon feedback bias in equation at 2 = 0, 0.5, 1 and 2. 
The vertical axis is the power spectra with baryons divided by power spectra 
with d ark matter only. The dots are simulation resultsofjvai^aalei^^l 
j2^n|) and the solid curves are the fitting formula of lHamois-Derans et al, 

<20l4 . 


4.2 Including Baryonic Effects 

In this subsection, we discuss the baryon-cooling effects on the 
matter power spectrum and the lensing dispersion. In our anal¬ 
ysis, we adopt previous results of high-resolution hydrodynamic 
simulations to incorporate the baryonic effects into the lensing dis¬ 
persion calculation. The baryon cooling naively enhances the mat¬ 
ter power spectrum by some tens of percent at fc = lO/iMpc”^ 
and the enhancement is more significant for smaller scales (e.g. 
Ijing et ^l20Q6l : iRudd etT]|200a : ICasarini et alJl2oTlh . Recently, 


two gr oups of the Overwhelmingly Large Simu lations (OWLS) 
project I Scha^e_et alj2010: |van Daalen et al.l20llh and the Illustris 
project I Vogelsberger et anbouJbl) suggested that the AGN feed¬ 
back suppressed Pm(fc) at intermediate scale (fc = 1 — 10/iMpc“^) 
but the baryon cooling strongly enhanced Pm{k) by orders of mag¬ 
nitude at small scales (fc > lOO/iMpc”^). These previous nu¬ 
merical results are roughly consistent but slightly different from 
one another, since it strongly depends on modelling of compli- 
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Figure 6. Matter power spectra with baryons (dotted curves) and without 
it (solid curves) aX z = 0 and 1. The solid cuiwes are our Halofit results, 
and the dotted curves are products of the baryon bias factor o times our 
Halofit results. 


cated astrophysical process such as star formation, radiative cool¬ 
ing, supernova and AGN feedback. In this paper, we use the re¬ 
sults of the most realistic model (“AGN_WMAP7 ”) in the OWLS 
project. They used extended version of Gadgets dSpringellllOOSll 
which used the treePM algorithm to calculate the gravitational evo¬ 
lution and smoothed particle hydrodynamics to follow the evolution 
of gas particles. They employed 512® CDM particles and equal 
number of baryonic particles in a cubic box of 100/iMpc~® on 
a side. Their Nyquist wavenumber is 16/iMpc“^ and their soft¬ 
ening length (rsoft = 7.8/i~^kpc) corresponds to a wavenumber 
27r/rsc,ft ~ 800/iMpc“^. They provided tabulated data of the mat¬ 
ter powe r spectra up to fc ~ 500feMp c~® on their website 0. Very 
recentlv. lHamois-Derans et aP d2014h made a useful fitting formula 
of the matter power spectra with baryons obtained by the OWLS 
project. They fitted the ratio of the power spectra with baryons 
^m,bary(fc, z) to that with dark matter only Pm,dmo(fc, z), called 
“the baryon feedback bias”, which is defined as. 


bl^ry{k,z) 


Tnijbary (li-, 

Tin,dmo(L, z) 


(13) 


FigS plots the baryon feedback bias in Eq. (H as a function of 
fc at z = 0, 0.5,1 and 2. As clearly seen in the figure, the AGN 
feedback suppresses the matter power spectrum at fc « 10/iMpc“® 
at some tens of percent, but the baryon cooling strongly enhances it 
at fc > 100/iMpc~® by an order of magnitude. The fitting formula 
can reproduce the numerical results within 5% relative error for 
k < 100/iMpc“^ at 2 < 1.5 dHamois-Derans et al ][20li) . 

In this paper, we simply multiply the bias factor (H to the re¬ 
sults Pm(fc) of our dark matter simulations to include the baryonic 
effects. The bias factor was obtained for WMAP 7yr cosmologica l 
model without the running spectral index l lvan Daalen et alJbol ih . 
However, we assume that the baryon feedback bias does not de¬ 
pend on the running parameters a, /3. Positive running a, /3 model 
probably enhances formations of small-scale clumps (such as dwarf 
galaxies and star clusters) at high redshift, which would affect the 
non-linear power spectrum. However, the dependence of the bias 


® http://vdll.strw.leidenuniv.nl/ 


on a, /3 is unknown at present and further investigations (hydro- 
dynamical simulations with running spectral index) are necessary. 
Fig|6] shows the power spectra with baryons (dotted curves) and 
without it (solid curves) at z = 0 (lower panel) and 2 = 1 (upper 
one). In this figure, Pm,bary(fc)fc^ is a monotonically increasing 
function for large k{> 100/iMpc“®)Q. Thus, the integration in the 
lensing dispersion ® strongly depends on fcmax (it does not for¬ 
mally converge at fcmax —i oo). In what follows, we will present 
the lensing dispersion results for several fcmax- 

Fig|7] shows the lensing dispersions with the baryons. Three 
panels are for various maximum wavenumbers in the integration 
fcmax = 300, 500 and 1000/iMpc“® (where the maximum 
fcmax is set to be 1000/iMpc“^, because the fitting formula of the 
bias factor bhar y becomes negative for k > 1000/iMpc“® at high 
redshift 2 > 2 dHamois-Derans et al.ll2014h ). As expected, the re¬ 
sults strongly depend on fcmax- Depending on the chosen kmax var¬ 
ious values of a, j3 are excluded, by the conservative (plausible) 
upper bound Gm^z ^ 1) ^ 0.12 (0.095). 

FiglU is the contour plot of the lensing dispersion with the 
baryons at 2 a = 1. The result Gm for fcmax = 300/iMpc“® 
is similar to the previous result in Fig|4l but Cm is much larger 
than the previous results for fcmax = 500 and 1000 /iMpc“^. The 
SNLS l-(2-)sigma constraints can be rewritten as a + 1.42)1 < 
0.145(0.222), 0.075(0.139) and 0.0006(0.054) for fc^ax = 
300, 500 and 1000/iMpc“®, respectively. Our numerical results of 

Gm{Zs\a,P,kma.fl at 0.6 ^ 2a < 2 and lOO/lMpC”^ < fcmax < 

1000/iMpc“^ can be fitted as. 


where x = 3.32a -|- 4.72/3 and 


mag 

(14) 


/o(Zs. fcmax) = [0.051 -f 0.009(2a - 1) - 0.009(2a - 1)^] 
X (1 -I- 0.09fcmax,lOo) , 

/l(2s, fc,nax) = [0.046 -b 0.017(2a - 1)] 

X (l -b 5.38fcmax,100 — 0.27fcnjax,10o) i 

(15) 


where fcmax,ioo = (fcmax/100/iMpc“^). The fitting function 
agrees with our numerical results within a relative error of 8.7% 
for 0 ^ a, 0.1. A byproduct of this analysis, is an easy deriva¬ 
tion of the validity of the baryon feedback. Taking a = /3 = 0, we 
hit the dispersion upper bound of g{z ^ 1) ^ 0.12 if we integrate 
up to fcmax = 1500/iMpc“®, wMle the ‘plausible’ cr < 0.095 is 
reached for fcmax = 950/iMpc“^. Hence the baryon bias has to 
be significantly modified and/or non-Gaussian effects have to be 
accounted for to explain the absence of detection of SNe lensing. 


5 CONCLUSIONS AND DISCUSSION 

We discussed the observational constraints on the small-scale 
power spectrum of the primordial curvature perturbation from SNe 
lensing dispersion. We modelled the common initial power spec¬ 
trum with running parameters a and j3 of the spectral index, and ran 


^ The original result of Pm,ba,iy{k, z)k^ in the OWLS project is not an 
increasing function of k, rather roughly constant at k > 30/rMpc-L Their 
dark-matter simulation result dmo(^)^) is slightly smaller than ours 
at small scales {k > lOO/iMpc”^), which results in the difference of 

^m,bary(^) ^)- 
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Figure 7. Same as Fig[2 but including baryonic effects on the matter power spectrum. Each panel corresponds to various maximum wave numbers in the 
integral of the lensing dispersion 0, fcmax = 300, 500 and 1000 hMpc“^ from left to right. The baryons clearly enhance the lensing dispersion compared 
to the previous results in Fig[3] 



Figure 8. Same as Fig|4] but including baryonic effects on the matter power spectra. Each panel is for fcmax = 300, 500 and 1000 hM.pc~^ from left to right. 
Note again that the thick solid line s correspond to 1- and 2-sigma observational upper bounds from JIO and JLA, while the white contours are the PLANCK 
jPlanck Collaboration et alj2014bh 1 and 2 sigma likelihood contours. 


the cosmological high-resolution N-body simulations to follow its 
non-linear evolution. We found that the small-scale enhancement 
in the primordial power spectrum is largely erased at late times 
during non-linear evolution, which suggests that it would be diffi¬ 
cult to probe the small-scale primordial fluctuations using late time 
measurements. For the dark matter models, the observational con¬ 
straints on a and P from the current observational data of the SNLS 
of JLA are weaker than the PLANCK results. However, if we in¬ 
clude the baryonic effects, which strongly enhance the non-linear 
matter power spectrum at small-scales (fc > 100/iMpc), the con¬ 
straints are comparable to or tighter than PLANCK. 

Let us discuss the accuracy of our Halofit formula in the high- 
k limit. The fitting formula was made to ht the simulation results in 
a finite range of 2/iMpc“^ < k < 600/iMpc“^ (see Appendix A). 
In the high-fc limit, the simulation resolution limits the accuracy. In 
the initial particle distribution, the linear Pm(fc) is included up to 
the Nyquist wavenumber IcNyq and is exactly zero for k > IcNyq, 
and thus our simulation results may underestimate the non-linear 


Pm(fc) for k > ^Nyq at late times. Evidence for such an under¬ 
estimation is the enhanced power spectrum in the high resolution 
simulations compared to the low resolution ones (IcNyq = 64 and 
40/iMpc“^, respectively). Additional hint for such an underes¬ 
timation is that at some point the nonlinear power spectrum be¬ 
comes smaller than the linear one for very large wavenumbers. 
For example at fc ~ 1000/iMpc“^ for z = 0, a = — 0.03 

and at fc ~ 120/iMpc“^ for z = 0, a = P = 0.06. This 
lack of resolution would degrade the accuracy at small-scales es¬ 
pecially for large a, P models. However, as seen in Fig[2l the high- 
and low-resolution simulations with aforementioned Nyquist wave 
numbers, give similar results with less than 20% difference for 
k < 300/iMpc“^, suggesting the safe convergence of our sim¬ 
ulations. Further studies investigating the non-linear Pm(fc) with 
positive Q, /3 at such smaller-scales are necessary. 

As we have seen, baryonic effects strongly influence our con¬ 
clusions. Baryonic simulation results strongly depend on assump¬ 
tions about various astrophysical processes (star formation and 
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AGN feedback, stellar mass function and so on ). Our results are 
based on two assumptions: i) We have used the Ivan Daalen et ^ 
(l201ib results. (In other words, we assume that the simulation of 
Ivan Daalen et alj ( 1201 ih is correct.) ii) The baryon feedback bias 
in (ES can be used even with large a and /3. Then, under these as¬ 
sumptions, we can safely set the strong constraints. It may very well 
be that assumption ii) is incorrect, and simulations of baryons with 
initial a and /? will produce a better behaved spectrum, similar to 
our dark matter simulation results. If our assumptions will be veri¬ 
fied, it is clear that baryons have an overwhelming effect, as they in¬ 
duce 0{1) changes of the power spectrum at fc > 5 /iMpc“^. Fur¬ 
thermore, the changes in the power spectrum are non-monotonic in 
fc. Last, with current simulation data, baryon cooling destroys the 
UV convergence, Pm{k)k^ diverges at fc ^ oo, making the lens¬ 
ing dispersion a UV sensitive quantity. It is therefore senseless to 
derive conclusions on the power spectrum without the inclusion of 
baryons. Moreover, any (late-time) astrophysical probe that is sen¬ 
sitive to such scales, requires an accurate baryon modelling before 
deriving any conclusions about cosmology. As discussed in sec l4.2l 
we need high-resolution hydrodynamic simulations with various 
initial conditions, that can probe very small-scales to reach a firm 
conclusion. Given the huge effect of baryons on lensing dispersion, 
SNe lensing might be a useful calibrator for baryonic processes. 
Especially since inclusion of baryonic effects in the Concordance 
Model {a = ft — 0), implies a diverging power spectrum, that is 
unphysical, and gives a lensing dispersion very close to the obser¬ 
vational bound for kmax ~ 1000, /iMpc“^. 

Another natural development is considering higher moments 
of the lensing pdf. Including small-scale fluctuations, the prob¬ 
ability distribution of the magnification would be strongly non- 
Gaussian (highly skewed). In fact, according to previous ray- 
tracing simulations in standa rd ACDM universe (e.g. lHilbert et alj 
l2007LTihkahashi et al.ll201ll) . the magnification PDF is rather a 
(modified) lognormal distribution with power-law tail at high /j.. 
Hence, the simple dispersion would not be sufficient to describe the 
lensing effects, and then the higher-order statistics such as skew¬ 
ness and kurtosis provide independent useful information. 

Finally, let us comment on future constraints on the running 
parameters a and ft, and on the power spectrum at small scales. 
Current tightest constraint on a and ft comes from PLANCK, 
which constrains them within 0(0.01) accuracy. Adding to the 
CMB observations, future radio telescopes measuring 21cm fluc¬ 
tuations such as the Square Kilometer Arra5l3 and the O mniscopj^ 
will constrain them within 0(0.001) accuracy (e.g. iMao et al 
l2008l : lKohri et al.l2013l : IShimabukuro et al.l2014l) . Small-scale per¬ 
turbations generate spectral distortions of CMB by Silk dumping 
due to photon diffusion, and future planned experiment PIXIE are 
also expected to probe the power spectrum up to sc ales as small 
as fc = lO'^Mpc"^, and potentially mea sure Q.,ft. dChluba et alj 
l2012l : |Powellll20l3 : IChluba & Jeongl2014l) . Given the baryonic ef¬ 
fects that swamp late time measurements, it seems that such early 
universe measurements are more promising in that respect. 
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APPENDIX A: FITTING FORMULA FOR NON-LINEAR 
MATTER POWER SPECTRA 


dex. Our formula is b ased on the halofit model dSmith et al.ll2003l : 
iTakahashi et al.ll20l^ hereafter T12), but slightly modified for the 
running models. 

To find the best fitting parameters in the halofit model, we use 
the standard chi-square analysis, which is defined as 


2 

X = 


= EE E 

1 . Z—Ofe —fcmi 


z') Pi, 

Pi,z)^ 


(Al) 


where Pi,model is the power spectrum in the theoretical model. 
Pi,sim is those in the simulation results, and i denotes the seven 
models of the running parameters (a, f)). The -)d is summed over 
redshifts z = 0, 0.3, 0.6,1 and 1.5. We use the wavenumber k 
larger than 2 /iMpc“^(= fcmin) where the Gaussian error of P(fc) 
is less than 1%. We set the maximum wavenumber fcmax so that 
the measured power spectrum is much larger (10 times larger) 
than the shot noise. Thus, fcmax depends on the parameters of a, j3 
and 2 . The maximum (minimum) fcmax is 600(100)/iMpc“^ for 
a = P = 0.06(0) atz = 0(1.5). 

First, we fit the simulation results for the no running model 
(a = /3 = 0). Here, we slightly modify the formula in T12. The 
halofit in T12 can be used up to k < SOfiMpc ^ within 10% 
error, but it asymptotically overestimate the Pm(fc) at very large 
fc(> 100/iMpc“^). The fitting parameters are the same as in T12, 
except for the following two parameters: 

logio cn = -0.3370 -F 2.1842neff -F 0.9604nefr -F 2.3764(7, 

7 n = 0.1027 - 0.0219neff -F 2.6785(7 

- (0.2866 -F 0.1206neff) In ( - ■ (A2) 

V /iMpc / 

In the last term, if one uses the wavenumber in unit of Mpc“^ (not 
fiMpc”^), then one has to insert a specific value of h. The equation 
( lA2t gives a minor correction to the formula in T12, which corrects 
the overestimation in T12 at fc > lOO/iMpc”^. The Pm(fc) with 
and without the correction in equation | |A2I are different within 
17%forfc < 100/iMpc“^ 

The ratios of the power spectra with the running to that with¬ 
out it are fitted as. 


P{k,z\a,P) 
P{k, 2; a = /9 = 0) 

with 


1-F 0.123 


/ k 


y/iMpc 


p{z-,a,P) 


(A3) 


p(2; a, /3) = (3.32q! -F 4.72/3) (1 -F (A4) 

The RHS of equation JA3b corresponds to an enhancement factor. 
Note that the parameters rieff and (7 in equation dA2b are evalu¬ 
ated in the no running model even when computing Pmik) for 
the running models. Using our fitting formula, the simulation re¬ 
sults can be reproduced within a relative error of 18.3(28.3)% for 
k < 300(600)/iMpc“^ at 2 = 0 — 1.5. The root-mean-square de¬ 
viation between the theoretical model and the simulation results is 
about 5.1%. 

Fig lAll shows the ratio of the power spectrum with the run¬ 
ning, Pm(fc, z; a, /3), to that without it, Pm(fc, 2 ; a = 0, /3 = 0). 
The red solid curves show the enhancement of power spectra de¬ 
scribed by our fitting factor in Ea. jA3b . which reproduces our sim¬ 
ulation results very well. 


This Appendix presents our fitting formula for the matter power 
spectra with running and running of running of the spectral in- 
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Figure Al. Ratio of the power spectra with the running Pm(fc, 2 ; a, /3) to that without the running Pm(k, z] a = 0, ^ = 0). The blue diamonds (black dots) 
are the ratios in the high-(low-)resolution simulations. The dotted curves are the ratios in the linear power spectra. The red solid cuiwes are our enhancement 
factor in equation <A3) . The small dumping at very high k(> lOOfeMpc ^) of the simulation results are due to the shot noise in the denominator Pm {k,z\a = 
0,/3 = 0). 
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